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ABSTRACT 



Context. The supernova remnant Cassiopeia A is a prime candidate for accelerating cosmic ray protons and ions. Gamma rays have 

been observed at GeV and TeV energies, which indicates hadronic interactions, but they could also be caused by inverse-Compton 

scattering of low-energy photons by accelerated electrons. 

Aims. We seek to predict the flux of nuclear de-excitation lines from Cas A through lower-energy cosmic rays and to compare it with 

COMPTEL measurements. 

Methods. Assuming a hadronic origin of the high-energy emission, we extrapolate the cosmic ray spectrum down to energies of 

10 MeV, taking into account an equilibrium power-law momentum spectrum with a constant slope. We then calculate the nuclear line 

spectrum of Cassiopeia A, considering the most prominent chemical elements in the MeV band and their abundances as determined 

by X-ray spectroscopy. 

Results. We show that the predicted line spectrum is close to the level of the COMPTEL sensitivity and agrees with conservative 

upper limits. 
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1. Introduction 

Cassiopeia A (Cas A) is one of the youngest known Galactic 
supernova remnants (SNR). Extensive observations in the radio, 
infrared, optical, and X-ray bands lead to an estimate d explosion 
date o f around 1680 AD (.Thorstensen et al...2001; jFesen et al.l 
l2006l) . Being the brighte st radio source in o ur Galaxy, it was 
first discovered in 1948 ("R vle & Srniflil ll948'). Its distance can 
be estimated to 3.4^^ ] kpc based on the combination of Doppler 
shifts and proper motions, whereas the angula r size of 2.5 ' cor- 
responds to a physical size of 2.34 pc (.Baars et al.lll977h . The 
optical spectrum obtained from distant Ught echos of the orig- 
inal blast indicates that Cas A was a type lib supernova of a 
~ 15 Mq main-sequence star and originated from the collapse 
of the helium core of a red supergia nt that had lost mos t of its 
hydrogen envelope before exploding dKrause et al.ll2008h . 

The blast wave drives an initially outgoing forward shock 
that can be seen as a thin X-ray edge expanding at roughly 
5000 km s"' . The interaction with circumstellar material and the 
interstellar medium results in a reverse shock that is driven back 
into the outgoing ejecta and expanding at roughly half the rate 
of the forward shock (Gotthelf et al. 2001; DeLanev & Rudnick 



12003 ). Emission at most wavelengths is dominated by a 30" 
thick "bright ring" where the ejecta are heated and ionized 
when they encounter Cas A's reverse shock. Viewed in X- 
rays, this shell consists of undiluted ejecta rich in O, Si, 
S, and Fe (Willingale et al. 2002; Laming & Hwang 120031; 
iLazendic et al . 2006). Faint X-ray filaments outside of the shell 
mark the location of the forward shock where nonthermal X- 
ray sy nchroton radiation is produced by shock-acc elerated elec- 
trons (lGotthelfetal."200l'; Vink & Laming"2003'). Detected by 
HEGRA (Aharonianetal. 2001), MAGIC (Albert etal. 200J) 
and VERITAS (Humensky 20oi), Cas A was the first SNR 
verified in TeV gamma rays. Recent observations with Fermi- 



LAT in the GeV range do not rule out either a leptonic or 
a hadronic emission scenario: A combination of nonthermal 
bremsstrahlung and inverse-Compton emission on the one hand 
as well as neutral pion decays on the other hand cou ld be respon- 
sible for the measured spectrum dAbdo et al.ll2010h . 

It has long been suggested that the shock waves associ- 
ated with supernovae could be sites of acceleration of cosmic- 
ray particles (Baade & Zwickv 1934). The combination of diffu- 
sive shock acceleration (DSA) processes in SNR shocks with 
transport effects in our galaxy can in theory produce the ob- 
served power- law spectrum of cosmic rays up to the "knee" of 
about lO'^ eV (Krvmskii"l97r;'Axford et al.'l977";'Bell"l978a"^ 



Bla ndford & Ostrik er 1978; Drury 1983; Jones & ElHson 1991 
To be consistent with the energy density of Galactic cosmic rays, 
the transfer of kinetic energy released in supernova explosions 
to cosmic rays must take plac e with an efficiency of ~ 10% 
(iGinzburg & Svrovatski i 1969). According to the DSA model, 
the charged particles scatter off irregularities in the magnetic 
field and increase their momentum by a fraction of the shock 
velocity v/c each round trip from the downstream to the up- 
stream region and back again. Because of the feedback of the 
accelerated particles on the spatial profile of the flow velocity 
and therefore on the particle distribution itself, the DSA pro- 
cess can be significantly nonlinear High-energy particles with 
considerable diffusion lengths sample a broader portion of the 
flow velocity profile and therefore experience a greater change in 
compression ratio than low-energy particles. That is why the par- 
ticle spectrum and the resulting photon spe ctrum tend to flatten 
with increasing energy (Baring et al. 19991; lEllison et al] 120001: 
iBell & Luce feOOlLlEliison et al., 2005). The correspondence of 
electron and proton spectra at high energies (Ellison et al. 200y) 
renders it extremely difficult to prove the acceleration of cosmic 
ray protons at SNR shocks directly. Though the presence of cos- 
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mic ray protons could be inferred by gamma rays resulting from 
collisions with ambient gas and subsequent pion decays, TeV 
gamma rays can also be produced by inverse-Compton scatter- 
ing of cosmic microwave background photons with the acceler- 
ated cosmic -ray electrons in SNRs. 

The detection of nuclear de-excitation gamma-ray lines 
could point a way out of this dilemma: The line spectra pro- 
duced in energetic collisions of cosmic rays with the ambient 
medium provide an unambiguous proof for the hadronic sce- 
nario of cosmic ray acceleration in SNRs. In contrast to other 
proposed scenarios (cf. Bozhokin & Bvkov 1997), we expect the 
excitation processes to happen near the proton acceleration site, 
so effects through Coulomb losses regarding the proton spectrum 
can be neglected here. The following sections are devoted to a 
short summary concerning the theoretical concepts of nuclear 
de-excitation processes as well as the derivation of the resulting 
gamma-ray fluxes in the MeV range with respect to current and 
future detection limits. Finally, we discuss the results and their 
implications in the last section. 



2. Theoretical framework 

Inelastic scattering of energetic particles on heavier nuclei as 
well as spallation reactions leaving the product nucleus in an 
excited state are followed by the emission of gamma rays in the 
range from 1 to 20 MeV through de-excitation processes. The 
spectral structure of these gamma rays is determined both by 
the composition and the energy spectrum of the energetic par- 
ticles and by the respective properties of the ambient medium. 
Especially cosmic rays with energies less than 100 MeV are 
suited well to be studied in gamma rays, because for cosmic 
rays with greater energies the gamma fluxes owing to p-p and 
p-a collisions followed by n^ desintegration are expected to 
be higher t han the gamma fluxes re sulting from nuclear de- 
excitations ("Men eguzzi & Reeveslll97 5). That is why the obser- 
vation of gamma-ray lines below 1 00 MeV offers the opportunity 
of studying astrophysical processes in great detail and of reveal- 
ing the origin of hadronic cosmic rays in SNRs. 

The basic ingredients in determining the profile of a gamma- 
ray line from energetic particle interactions can be summarized 
as follows: According to Ramatv et al. ( 1979), the probability of 
photon emission per second into solid angle d{cos 9q) d<po can be 
written as (the z-axis is chosen in the direction of the incident 
particle) 



and E from the equations 



dcr 
dPy = niV—-{E, el) d(cos 9*.) dcpr x 

X g{E, e*, Oo, 4>r - </>()) <i(cos 6*0) #0- 



(1) 



In the center-of-mass frame, the interaction produces an excited 
nucleus with recoil velocity in ii(cos 0*)d(f>r- E represents the en- 
ergy of the incident particle before the interaction, «, the number 
density of the target particles, v is the incident particle's velocity, 
dcr/dO.* is the center-of-mass differential cross section, and g is 
the angular distribution of the resulting gamma rays, (pr and 0o 
are azimuth angles measured in the (x, 3')-plane, 0* and 6*0 are po- 
lar angles given with respect to the z-axis. To finally determine 
the gamma-ray spectrum, Eq. ^ can be integrated over cosO*, 
(f>f, cos 0() and E by using the Monte Carlo simulation technique. 
Choosing random numbers (uniformly distributed from to 1), 
the integrations can be carried out by solving for cos 0*, 0,., cos Oq 



o-(E) J dU.* 



Ri^—, 

In 



R^^C I vNp{E')o-{E')dE', 



Ra^ -(l+cos^o)- 



(2) 
(3) 
(4) 
(5) 



Here Np(E) is the number of incident particles per unit energy, C 
is a normalization constant. The probability of observing gamma 
rays of energies between Ey and Ey + AEy is then proportional 
to the sum of all angular distributions g for which Ey is in range. 
For a detailed description of the outlined methods above and 
a deeper insight into the different reaction types as well as the 
de rivation of lin e produ ction cross sections, we refe r the reader 
to 'Rama tv et aP (1 19791) and iKozlovskv et al.1 (l2002h and refer- 
ences therein. 



3. Results 

To validate the chances for success of the introduced approach 
concerning Cas A, we first study the exemplary case of the 
carbon line at 4.4 MeV. Considering a scenario in which the 
gamma-ray emission of Cas A is modeled by a hadronic fit based 
on the 7r"-decay of accelerated hadrons, the best-fit proton accel- 
eration spectrum is given by Qp(p) o^ p^^'^. The resulting pro- 
ton energy content of H'p = J^Q^^^^^MpPdp = 4 X 10'*'' erg 
corresponds to ~ 2 % of the estimated SNR kinetic energy 
(lAbdo et al.l l2010h . By extrapolating the high-energy proton 
spectrum down to the MeV-range, the gamma-ray flux emitted 
in the 4.4 MeV nuclear de-excitation line can be approximately 
calculated from 



1 



471^2 



/ 



"c Qp(p)o-(p)v{p)dp, 



(6) 



where rtc ~ 10 cm -^ is the adopte d mean density of carbo n 
atoms in the interaction region (cf. iLaming & Hwan^ l2003h . 
V the velocity of the accelerated protons and cr the cross sec- 
tion for the inelastic scattering processes. At first glance, line- 
broadening effects or additional contributions from unresolved 
gamma-ray lines in heavy nuclei and lines from long-term ra- 
dioactive nuclei are neglected here. Using d - 3.4 kpc (see 
above) and t he cross section for the reaction ^^C{p,p'y^C* 
given by Ra matv et al.l d 19791) yields a flux of ~ lO^^cm ^s ' 
at 4.4 MeV. This flux value is close to the sensitivity limi t of the 
COMP TEL experiment: Following the analysis of Strong et al.l 
(l2000h concerning the MeV continuum emission from Cas A, 
only an upper limit of 1.4 x lO^^cm^^s"' was obtained in the 
3 - 10 MeV energy range. According to Ivudin et al. (1995), the 
COMPTEL line sensitivity is indicated with ~ 10"^ cm^^s"'. So 
only a future gamma-ray mission with enhanced sensitivity in 
the MeV range will be able to obtain final results concerning the 
detection of de-excitation lines in Cas A. 

To compute the whole nuclear de-excitation spectrum for 
the specific case of Cas A , we u sed the Monte-Carlo code de- 
veloped by iRamatv et al.1 (Il979ln (see also Sect. |2]i. Besides 



ihttp : //Iheawww . gsf c . nasa . gov/users/ramaty/code . htm| 
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Fig. 1. Calculated gamma-ray spectrum for the specific case of Cas A using the assumptions described in the text. 10^ photons are 
binned into energy intervals of widths ranging from 2 to 5 keV as described in Ramaty et al. ( 1979). For example, the '^C line at 
4.4 MeV has a flux that is comparable to the sensitivity of COMPTEL. Note that the continuum component caused by nonthermal 
electron bremsstrahlung is not taken into consideration here. 



the ingredients already mentioned in the previous paragraph, 
the calculations take into account the following assumptions: 
The acceleration scenario for cosmic rays is assigned to the 
reverse-shock side and the chemical composition of th e accel - 
erated cosmic rays is inferred from lEngelmann et alj (Il990h . 
The composition of the ambient gas, in fact a mixture of both 
massive Wolf-Rayet winds and subsequent supernova ejecta (cf. 
iLingenfelter & Hi gdon 2007), is described by the use of re- 
sults from X-ray spectroscopy (cf. Table [T]). The abundances 
of H and He are deduced from optical measurements by 
IChevalier & Kirshnerl ( Il979h . The resulting total mass is in line 
with the Wolf-Rayet-supernova scenario, i.e. there is no room for 
additional amounts of hydrogen that would enhance the pion vs. 
the nuclear de-excitation yields. Unresolved gamma-rays from 
heavy nuclei and lines from long-term radioactive nuclei are 
also included. Together with the consideration of recoiling tar- 
get particles, this leads to a significant broadening of the lines. 
Although some nonthermal emission is associated with the for- 
ward shock, recent studies showed that electron acceleration to 
multi-TeV energies is likely to take place mainly at the reverse 
shock within the supernova ejecta (Helder & Vink 2008), mak- 
ing Cas A to a unique object for studying particle acceleration 
at the reverse-shock side. The resulting spectrum is depicted in 
Fig.[T]and agrees well with the approximate calculation above. 



4. Discussion 

As can be seen in Fig. [1] the adopted acceleration scenario for 
cosmic rays in Cas A leads to a flux of nuclear de-excitation 



Table 1. Mean measured abundance mass ratios and rms scat- 



ter resp. upper limits according to the results of Willingale et aL 
(2002), Docenko & Sunyaey (2010) and iChevalier & Kirshnei 
(119791) . 



ratio 



H/Si 


< 2.29 X 10-5 


- 


He/Si 


< 4.93 X 10-3 


- 


C/Si 


1.76 


0.88 


O/Si 


1.69 


1.37 


Ne/Si 


0.24 


0.37 


Mg/Si 


0.16 


0.15 


S/Si 


1.25 


0.24 


Ar/Si 


1.38 


0.48 


Ca/Si 


1.46 


0.68 


FeL/Si 


0.19 


0.65 


FeK/Si 


0.60 


0.51 


Ni/Si 


1.67 


5.52 



Notes. The mass ratios are given relative to solar values. To faciliate the 
comparison, the data from IChevalier & Kirshnej dl979l) for H and He 
are indicated relative to Si, too. 



lines that would be clearly detectable by a gamma-ray telescope 
with enhanced sensitivity as sucessor t o the COMPTEL m ission 
(e.g. the proposed GRIPS mission bv iGreiner et al.! !2009). The 
line emissivity is additionally boosted because the thermal target 
gas reflects the heavily enriched abundances of the Wolf-Rayet 
progenitor star. Though the detailed line characteristics always 
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depend on the precise knowledge of the supernova ejecta's com- 
position, the natural process of element synthesis in the pro- 
genitor star and the supernova explosion mechanisms lead to 
peculiar properties in the gamma-ray spectrum. A unique fea- 
ture arises because the C and O lines in the 4-6 MeV band 
dominate the line flux from the Ne-Fe group in the 1-3 MeV 
band. Of course, the reasonability of extrapolating the proton 
spectrum obtained from high-energy measurements down to the 
MeV range has to be discussed: It is commonly believed that 
strong Coulomb losses would quench the particle spectra of non- 
thermal particles below GeV energies, but the steep spectra of 
nonthermal particles in solar flares may tell a different story. 
Furthermore, the results of modeling the acceleration process 
of cosmic rays and the evolution of SNRs sh ould be taken into 
account: According to the kinetic models of i Berezhko & VolM 
(Il997h andHerezhko & Volk (2000), the overall proton momen- 
tum spectrum at the shock front is close to a pure power law 
in the entire momentum range down to MeV at low injection 
rates, for high injection rates the low-energy part {p « mc) 
of the overall cosmic ray spectrum is slightly steeper than the 
high-energy part. The application of the kinetic model approach 
with respect to the spe cific case of Cas A confirms these results 
(IBerezhko et al.ll2003h . So the assumption of an in-situ cosmic 
ray spectrum that is as steep as at higher energies and a localiza- 
tion of the inelastic scattering processes near the reverse shock 
site seem to be quite reasonable and theoretically motivated. 

5. Conclusion 

We have shown that if cosmic rays are accelerated at the reverse 
shock in the Wolf-Ray et supernova remnant Cas A, as indicated 
by its strong GeV and TeV emission, the resulting flux of car- 
bon and oxygen de-excitation lines is marginally detectable with 
COMPTEL. Preci sion spectroscopy with future MeV-missions 
such as GRIPS ( Greineret"aLll2009i) will put the theory of cos- 
mic ray acceleration to a crucial test and will permit us to de- 
termine the overall efficiency of cosmic ray acceleration. As a 
corollary, these measurements would also determine the yields 
of spallation products from nuclear collisions that affect the 
abundances of the cosmologically relevant light elements such 
as Li, Be, and B. It must be noted that the prediction of the line 
flux rests on an extrapolation of the cosmic ray spectrum de- 
termined at high energies, which is in line with difiiisive shock 
acceleration theory, but which is still subject to major system- 
atic uncertainties associated to nonlinear effects such as shock 
broadening caused by back reactions. 
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